DYNAMICS AND THERMODYNAMICS OF PHASE TRANSITION IN 

HOT NUCLEI 
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Abstract 

The dynamics and thermodynamics of phase ttansition in hot nuclei are studied through experimental 
results on multifragmentation of heavy systems (A> 200) formed in centtal heavy ion collisions. Dif- 
ferent signals indicative of a phase transition studied in the INDRA collaboration are presented and their 
consistency is stressed. 
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1. INTRODUCTION 

Liquid-gas type phase transitions are com- 
monly observed in systems with short-range re- 
pulsive and longer-range attractive forces, such 
as macroscopic fluids with the van der Waals in- 
teraction. The nuclear equation of state is very 
similar to that of non-ideal gases, which allows 
to foresee the existence of different phases of nu- 
clear matter. But can one define phase transi- 
tions in "small" objects, when the dimension is 
of the same order of magnitude as the range of 
the force which binds them? In this sense nu- 
clei as well as clusters of galaxies are small. Re- 
cent works state that statistical mechanics based 
on the Boltzmann' definition of entropy (S = 
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k. log PL) allows to define phase transitions in 
small systems [[]]]. The nucleus at zero or mod- 
erate temperature, because of its quantal nature, 
is assimilated to a liquid phase. A nuclear gas 
phase was characterized, for instance, through 
the properties of vaporised quasi-projectiles from 
95 AMeV Ar+Ni reactions, fully predicted by 
modelling a van der Waals gas of fermions and 
bosons in thermodynamic equilibrium [0]. The 
liquid-gas coexistence region was thus naturally 
connected to multifragmentation, break-up of a 
nuclear system in several fragments of various 
sizes. 

Most of the works on nucleus phase transition 
relied on static properties such as the charge (or 
mass) partitions measured in nuclear reactions. 
The observation of a power law on a limited 
range of fragment charge measured in inclusive 
experiments was claimed to be consistent with a 
critical behaviour, according to Fisher's droplet 
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model [JSp, and also to percolation models [Q]. 
Recent and well selected data from the ISIS col- 
laboration were analysed in this framework, and 
a liquid-gas coexistence curve was derived |J]. 
In the same line, critical exponents were reported 
to have measured values in agreement with those 
of a liquid-gas model [0, ||] . 

Caloric curves relate the energy and the temper- 
ature of nuclei. The observation of caloric curves 
showing a bending or a plateau, similar to what 
is observed for the boiling of water, was claimed 
as a proof of the nuclear phase transition [Q] . 

A new signal of phase transition of nuclei was 
found in the observation of a negative branch of 
the microcanonical heat capacity for Au quasi- 
projectiles excited between 3 and 6 AMeV JTC[]. 

This paper presents a compilation of the re- 
sults which deal with central collisions between 
heavy ions (A pro j + A tgt > 200), obtained so far 
in the INDRA collaboration. Most of them will 
be detailed in other contributions to this confer- 
ence. Different signals which may be connected 
with the occurrence of a phase transition in nuclei 
were found, each of them bringing a small stone 
to the construction of a consistent framework. 

2. PHASE TRANSITION: DYNAMICS AND 
THERMODYNAMICS 

2.1. Thermodynamics aspect 

As exemplified in fig. [1], where isotherms are 
drawn in a pressure-density plane, the EOS of 
nuclear matter resembles that of a van der Waals 
gas. A region of negative incompressibility, bor- 
dered by the spinodal line, and located inside the 
coexistence region, appears in this diagram. This 
spinodal region is a zone of mechanical insta- 
bility, where an increase of pressure leads to an 
expansion of the system; on entering this zone 
the system will tend to recover thermodynami- 
cal equilibrium by separating into two coexisting 
phases, a low-density (gas) phase and a higher 
density (liquid) phase. Note that the time for the 
system to go from homogeneity to inhomogene- 
ity is finite, in other words a phase transition has a 
dynamics (see next subsection). The spinodal re- 
gion is also thermally unstable, as the heat capac- 
ity, c = dE/dT, is negative inside this zone. Ac- 




P/Pc 

Fig. 1 : Equation of state relating the pressure and the den- 
sity (normalised to the critical values) in nuclear matter. 
The curves represent isotherms. The dashed-dotted line is 
the coexistence line and the dotted line the spinodal line. 

cording to this formula, it seems easy to evidence 
a negative heat capacity, which should appear, in 
a microcanonical frame, as a backbending of the 
caloric curve. It is not so simple in the case of 
nuclei, as their temperature, energy, density can 
only be modified by a more or less violent nuclear 
collision; one is dealing with isolated systems 
without external constraints, which eventually in- 
duce a phase transition in the course of their evo- 
lution in a (temperature, pressure, energy) space. 
The followed path is not known and there is no 
a priori reason why the resulting caloric curve 
would have a plateau, or a backbending jrT|] . 
Another and more robust signal was proposed 
to recognise a phase transition, abnormally large 
fluctuations of the kinetic part, E k of the energy 
(E = E k + E pot ) [jrj]. In this case the heat capac- 
ity can be expressed as c tot = c\/ (c k — Aal/T 2 ), 
with c k = d(E k )/dT. When the normalised fluc- 
tuations exceed the canonical value c k , the heat 
capacity becomes negative, as depicted in fig. 0. 
These variables may be determined in an experi- 
ment, provided that one can reconstruct the parti- 
tion present at the "freeze-out" stage, defined as 
the instant when nuclear interaction between the 
present species vanishes. 

Finally this static statistical physics approach 
supposes that the ensemble of collisions consid- 
ered constitutes a statistical ensemble. This en- 
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Fig. 2: Variation versus energy of the kinetic fluctuations 
and heat capacity (top) and of the total heat capacity (bot- 
tom). 

semble is characterized by global variables re- 
flecting the pertinent information resulting from 
the dynamics of the collisions or (and) the sorting 
of the events. Statistical physics concepts can be 
safely used providing that there is adequacy be- 
tween the ensemble treated and the applied for- 



malism [13] 



2.2. Dynamics of the phase transition 

Is it possible to explore the fugitive dynamics 
of the nuclear phase transition? A positive an- 
swer was recently given in the INDRA collab- 
oration. Heavy nuclear systems were produced 
at the same energy (~ 7 AMeV), but with very 
different total masses (393 and 248), through 
central collisions between 36 AMeV Gd and U 
and 32 AMeV Xe and Sn; the resulting mea- 
sured charge distributions superimposed while 
the fragment multiplicities scaled as the total sys- 
tem masses Ml 411. While this indicates the domi- 



nance of phase space effects, with stochastic frag- 
ment production governed by a Boltzmann fac- 
tor, it may also sign the presence of bulk insta- 
bilities inducing density fluctuations, and finally 
leading to the spinodal decomposition of the nu- 
clear system. And indeed spinodal decomposi- 
tion is, together with nucleation, one of the pro- 
cesses describing the dynamics of the liquid-gas 
phase transition in nuclear matter jJISQ. The most 
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Fig. 3: Density fluctuations and fragment formation in the 
spinodal region [|l6||. 

unstable mode(s) in the system induce homo- 
geneous density fluctuations which will be am- 
plified by the mean field with a characteristic 
time which decreases with the temperature of the 
system. Locally zones close to normal density 
may appear, leading to the formation of frag- 
ments whose radii and relative distances depend 
on the wavelength of the mode, as schematized 
in fig. [3[ In this scenario, one expects a narrow 
fragment charge distribution, loosely dependent 
on the total charge of the system provided it is 
large enough. In finite size hot nuclear systems 
several modes are equally probable, may beat, all 
effects which blur the initial picture and lead to 
much broader charge distributions. This is ex- 



32AMeVXe+Sn 
36AMeVGd+U 




Fig. 4: Charge distributions from central heavy ion colli- 
sions. The lines show the results of stochastic mean field 
simulations while the symbols display the experimental re- 
sults. Adapted from Jl^]. 

emplified in fig. ffl, where the lines display the 
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charge distributions obtained through a complete 
stochastic mean field simulation of central nu- 
clear collisions between heavy ions, the Brow- 
nian one-body dynamics (Bob): the mean field 
is complemented by a Brownian force in the ki- 
netic equations, which mimic a stochastic colli- 
sion term. For the two reactions considered here, 
36 AMeV Gd+ U, and 32 AMeV Xe+Sn, the 
system is driven into the spinodal region, where 
it multifragments through spinodal decomposi- 



tribution [|2"2"|], in order to make its statistical er- 



tion [ |18| ] . The de-excitation of the hot fragments 
so formed is followed simultaneously with their 
propagation in their Coulomb field. Finally the 
simulated events (Bob events) are passed through 
a software replica of the INDRA detector. Obvi- 
ously the final distributions (lines in fig. |]) are 
far from narrow, being continuously decreasing 
with increasing Z without any visible peak. But 
they perfectly reproduce the experimental dis- 
tributions measured in the two systems quoted 
above (symbols in fig. |]), particularly the inde- 
pendence of the charge distribution against the 
mass of the system is recovered [[T7[]. It is worth- 
while to mention that kinetic properties of the 
outgoing channel are also well accounted for, 
particularly the fragment-fragment velocity cor- 
relations JT^]. 

These simulations show that the phase transi- 
tion, or the spinodal decomposition process can 
not be evidenced by a single variable like a 
charge distribution. It might however well be 
that a few events kept a memory of the pro- 
cess, through a preferred break-up in equal-sized 
fragments. The charge correlation function de- 
fined as: C = Y cor (Z av ,az)/Y uncor (Z av ,az) is 
very efficient to enlighten any extra production 
of rare event types [|2"0"|], as it makes use of the 
full per-event information through two variables, 
the average charge Z av and the standard devia- 
tion az- Equal-sized fragment partitions (small 
oz), were looked for in Bob events from 32 
AMeV Xe+Sn central collisions. The variable 
Z av was replaced by Mf x Z av as events with 
different fragment multiplicities Mf were mixed 
to improve statistics [|21|]. The uncorrelated yield, 
Y U ncor(Z av , a z ), obtained in ref [|2(J |2l|] by mak- 
ing pseudo-events from fragments belonging to 
different events of the same sample, may also 
be analytically calculated from the charge dis- 



ror negligible everywhere. The correlation func- 
tion shown in fig. [5| (left) uses the analytical 
Yuncor an d the standard deviation of the mea- 
sure: o z = (££ - Z av f/Mffl\ at vari- 
ance with ref [|2(], |T[, The general trend is a de- 
crease of the correlation function with decreasing 
az, as expected from the broad charge distribu- 
tion, except in the smaller bin, az < 1, where it 
increases again. For 36 < M f x Z av < 60 this 
increase is statistically significant, with a con- 
fidence level larger than 98%. But only 0.6% 
of the events contribute to these enhanced par- 
titions, while all simulated events multifragment 
through spinodal decomposition. Due to this 
small number one may wonder about the real 
meaning of these peaks at small az, because their 
statistical significance must be weighted against 
a "background", corresponding to the extrapola- 
tion of the trend of the measured correlation func- 
tion for decreasing a z . This extrapolation was 
empirically performed with an exponential func- 
tion [ |2~T| ] . The method was recently improved by 
building the uncorrelated yield following an inde- 
pendent emission model under the constraint of 



charge conservation of the total system [[22]]. The 
result, on the same Bob event sample, is shown 
in the right part of fig. [5|. All correlations due to 
charge conservation are now suppressed, and the 
correlation function is equal to 1 wherever no ad- 
ditional physical correlation is present. To make 
the picture stronger, the peaks with a significance 
lower than 95% have been flattened out in fig. |5] 
(right). Again one observes peaks for small val- 
ues of az, the percentage of events concerned 
is 0.7%, which fully validates the published re- 
sults [glp. 



In short simulations of central heavy ion colli- 
sions show that fingerprints of the spinodal de- 
composition of a nuclear system can be found, 
namely a weakly but unambiguously enhanced 
proportion of equal-sized fragment partitions, 
and thus evidence its phase transition. A sta- 
tistical model, like any model just based on 
phase space, will not show this type of phe- 



nomenon [ 22] . 
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Fig. 5: Charge correlation function for events generated by stochastic mean field simulations of central 32 AMeV Xe+Sn 
collisions, without (left) and with (right) charge conservation constraint on the uncorrelated yields. Note the different 
orientation of the two figures. The charge axes are Mf x Z av on the left, and Z av on the right, while the uz axes are the 
same. In both cases events with 3 to 6 fragments are included. 



3. EXPERIMENTAL INDICATIONS OF 
MULTIFRAGMENTATION IN THE SPINODAL 
REGION 



Two reactions were studied by the INDRA col- 
laboration which lead to systems with close total 
mass and charge (Z~105, A~250). The symmet- 
ric Xe+Sn system was studied between 25 and 
50 AMeV, and a class of events characterized as 
multifragmentation of "fused" systems was rec- 
ognized at and above 32 AMeV. Event selection 
was performed on completeness and event shape 
criteria, and the multifragmentation of deformed 
objects was not considered in these samples [ jlTl 
[2~3"| ]. The asymmetric Ni+Au system was stud- 
ied from 32 to 90 AMeV. Owing to the thresh- 
olds of INDRA, and to the smaller cm. velocity 
for this asymmetric entrance channel, the num- 
ber of complete events is very small. Methods 
selecting classes of events were developed, less 
stringent on completeness requirements. Two 
of these methods were used for isolating events 
samples of central collisions with close distribu- 
tions of global observables: the Principal Com- 
ponent Analysis method [ |2"4"| , [2~5| ]) for the sample 
on which the heat capacity was determined and 
the Discriminant Factorial Analysis method [|25J]) 
for the sample used for charge correlations. 



3.1. The symmetric entrance channel 

The heat capacity and the charge correlations 
were both examined in the same samples of 
Xe+Sn central collisions. Charge correlations are 
presented in fig. for the four incident ener- 
gies [j2"7|]. The correlation functions shown here 
are calculated without the charge conservation 
constraint, and the figure only presents them for 
the first bin oz £ [0 — 1] . The data corresponding 
to an enhanced production of equal-sized frag- 
ment partitions with a high confidence level are 
indicated by darker points. This enhancement is 
only observed at 32 and 39 AMeV, and corre- 
spond to 0.2 and 0.4% of the selected samples. 
These percentages are close to those observed in 
the Bob simulations, and may be interpreted in 
stating that the dynamics of multifragmentation 
occurs by spinodal decomposition in all the se- 
lected central collisions. 

The heat capacity was determined on recon- 
structed freeze-out configurations out of the mea- 
sured data [j2~8[l. The excitation energy scale 
was obtained by calorimetry, including all mea- 
sured fragments and twice the light charged par- 
ticles emitted between 60 and 120deg in the cm. 
Neutrons were estimated from mass conserva- 
tion. The radial expansion energy was deter- 
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Fig. 6: Charge correlation functions (without charge con- 
servation constraint) for central collision events between 
Xe and Sn. The lines indicate the extrapolated background. 
The darker points correspond to event partitions enhanced 
with a high significance. From [p7[] 
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Fig. 7: Heat capacity versus the excitation energy (cor- 
rected for the collective energy) measured in central col- 
lisions between Xe and Sn. From [ p8| ] 

mined from a comparison of the data with re- 
sults of a statistical multifragmentation model, 
and was subtracted from the excitation energy to 
give the energy distributions shown by the lines 
in fig. [7|. The potential energy was taken as the 
sum of the mass balance and of the Coulomb en- 
ergy estimated through the Wigner-Seitz approx- 
imation, with a freeze-out volume V = 3 x Vq. 



The influence of these assumptions on the ampli- 
tude of the fluctuations is thoroughly discussed 
in []Z9[] . The heat capacity, derived as explained 



in sect. 2 2.1, is plotted as grey zones (error bars) 
in fig. [7|. Negative branches are observed below 
E* =6-6.5 AMeV at the four incident energies. 
At 32 and 39 AMeV there is a significant overlap 
between the energy distribution and the negative 
heat capacity; this indicates that multifragmenta- 
tion occurred in the spinodal region, in agreement 
with the signal given by the charge correlations. 
At and above 45 AMeV, both the non observation 
of enhanced equal-size fragment partitions and 
the positive value of the heat capacity plead for 
the system multifragmenting in the coexistence 
region outside the spinodal zone. 

3.2. The asymmetric entrance channel 




Fig. 8: Charge correlation function (without charge con- 
servation constraint) versus the standard deviation az for a 
given value Z av = 12. The line indicates the linear back- 
ground. From i fiot i 

Enhanced partitions of equal- sized fragments 
were looked for in central Ni+Au collisions at 32 
and 52 AMeV []3~T]]. A statistically significant en- 
hancement only appeared at 52 AMeV, as shown 
in fig. |3 3.2| , where the correlation function is cal- 
culated without charge conservation constraint. 
Results with the new method can be found in 
the contribution of B. Guiot, and agree with the 
above statement. 

The heat capacity was calculated for incident 
energies of 32, 52, 74 and 90 AMeV. Precisions 
on the implied assumptions are detailed in the 
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contribution of O. Lopez to this conference, in 
which the results are also unveiled. The energy 
scale is here obtained from calorimetry, using 
products emitted between 60 and 120deg cm. 
whatever their charge, due to stronger contribu- 
tions of forward preequilibrium emissions with 
respect to the symmetric entrance channel. A 
negative branch of the heat capacity appears be- 
low 6-7 AMeV at the four incident energies. The 
data at 32 AMeV should be taken with caution, 
as the average completeness of the events is poor 
(65% of the total charge). Indeed because of the 
low value of the average excitation energy (3.5 
AMeV), one would expect to observe the first di- 
vergence of the heat capacity, while only the sec- 
ond one is seen, in the higher part of the energy 
distribution. At the three other incident energies, 
the average excitation energy varies little (5.5, 7, 
7 AMeV) and the heat capacity is negative in the 
first half of the distribution. 

For this asymmetric entrance channel, multi- 
fragmentation in the spinodal region is shown 
by dynamical and thermodynamical signals at 52 
AMeV only. For collisions at higher energy only 
the negative heat capacity was looked for, and 
found. At 32 AMeV there is no sign of spin- 
odal decomposition, and the data resemble more 
a fusion-evaporation process [ffifl; surprisingly 
only the negative part of the heat capacity was 
observed despite the low excitation energies cov- 
ered. 



4. NEW SIGNALS EVIDENCING A PHASE 
TRANSITION 

4. 1. Universal fluctuations of an order parameter 

The recently developed theory of universal fluc- 
tuations in finite systems provides methods to 
characterize critical and off-critical behaviours, 
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without any equilibrium assumption D33Q- In this 
framework, universal A scaling laws of one of 
the order parameters, p, should be observed: 

(p} A P(p) = <j>((p-p mp )/(p} A ) 
where p mp stands for the most probable value 

of p. A=l/2 corresponds to small fluctuations, 
at ~ (p), and thus to an ordered phase. Con- 
versely A=l occurs for the largest fluctuations 
nature provides, a"t ~ (p) 2 , in a disordered phase. 
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Fig. 9: Normalised first and second moments of the scaled 
distributions of the event heaviest fragments, for Xe+Sn 
central collisions at five incident energies. The lines A = 1 
and 1/2 are shown to guide the eye. Note the two logarith- 
mic scales. Adapted from [P4I. 

Scaling laws may be expected for an order pa- 
rameter p increasing with time, for instance the 
fragment multiplicity in a fragmentation process, 
or the largest fragment size in an aggregation 
process. The method was applied to the Xe+Sn 
central collision samples described in sect. 
extended to the measurements performed at 25 
AMeV [^2p . The fragment multiplicity distribu- 
tions at the five energies scale with a small vari- 
ance (A=l/2), while the largest fragment distri- 
butions follow a A=l/2 scaling at low energy (25 
and 32 AMeV) and a A=l scaling for the three 
highest energies; fig. |^ illustrates the variation 
of 

Zmax vs - Zmax- Following the authors of 
this study, this indicates the transition from an or- 
dered phase (evaporation?) to a disordered phase, 
the fragments being produced following an ag- 
gregation scenario. Note that spinodal decompo- 
sition well enters this type of scenario. 

4.2. Bimodality of an order parameter 
An extended definition of first order phase tran- 



sitions in finite systems was proposed in [|34|], 
based on topology anomalies of the event distri- 
bution in the space of observation. In this frame- 
work when a nuclear system is in the coexis- 
tence region, the probability distribution of an or- 
der parameter is bimodal, each component being 
a phase. The Xe+Sn central collisions samples 
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Fig. 10: Distributions of the differences of total charges 
bound in big (Z>13) and small (Z E [3 — 12]) fragments 
in events from Xe+Sn central collisions. From fE7fl 



were also analysed this way. Having in mind the 
liquid-gas coexistence region, one has to guess 
which fragments belong to each of the phases. 
Owing to the detailed shape of the charge dis- 



tributions Q27[], Z=12 was tentatively put as a 
limit between the two phases. The chosen sorting 
parameter, (]T Z z > 13 - ]T Z Ze[3 _ 12] )/ J2 Z z > 3 , 
may be connected with the density difference of 
the two phases (pl — Pg), which is an order pa- 
rameter JT3|]. The results are displayed in fig. [TO: 



none of the distribution is Gaussian, which would 
sign a pure phase; but the gas phase is dominant 
above 45 AMeV, while the liquid phase is more 
important at 32 AMeV. The two components are 
roughly equal at 39 AMeV. 

Similar results were obtained on the Ni+Au 
central collision samples selected with the PCA. 
A similar sorting parameter, making use of only 
the three largest fragments, was used. Its proba- 
bility distribution exhibits two nearly equal peaks 
at 32 AMeV, a big "gas" peak with a small shoul- 
der at 52 AMeV and a single phase picture at 90 



AMeV [25] 



5. SYNTHESIS OF THE OBSERVATIONS 
FOR SYSTEMS OF 250 NUCLEONS 

A synthesis of the signals which may sign a 
liquid-gas phase transition for nuclei of mass 
~250 is presented in table I. The thermal and ra- 
dial energy scales come from SMM simulations, 
with respective uncertainties of 1 and 0.5 AMeV, 



depending, for instance, on the event selections, 
the freeze-out volume chosen . . . [ E41 E8p. The 
first evidence coming from this table is that for 
the asymmetric entrance channel, multifragmen- 
tation of the fused system formed in central colli- 
sion is not associated to a radial expansion. This 
comes at variance of other results for similar 
systems in the same energy range, for instance 



Kr+Au [[35J. Note that the result comes directly 
out of the data, as for Xe+Sn at 50 AMeV and 
Ni+Au at 90 AMeV, the partitions are similar, 
while the average kinetic energies of a given frag- 
ment strongly differ, and are constant for Z be- 



tween 12 and 30 for the latter system p5| ] 

For both entrance channels, the systems reach 
the spinodal region and remain there for a time 
sufficient to allow the development of spinodal 
instabilities for thermal energies between 5 and 
6.5 AMeV, as attested by the observation of 
favoured equal-sized fragment partitions and the 
negative values of the heat capacities. When the 
thermal energy increases, the expansion energy 
seems to have a decisive role: without expansion, 
the picture is similar to the one just described, up 
to e t h =7.5 AMeV, while an expansion energy 
above 25% of the thermal energy would make 
the system cross too fast the spinodal region, to- 
ward the coexistence region or even the gas re- 
gion. In this case neither the spinodal decompo- 
sition nor the negative heat capacity are observed. 
The critical excitation energy, given by the max- 
imum thermal energy of the systems described 
in this paper and shown to have fragmented in 
the spinodal region is at least 7.5 AMeV, much 
above the value given in 0]. It must however 
be noted that the "multifragmentation" partitions 
analysed in this reference comprise essentially 
one large cluster and one intermediate mass frag- 
ment, together with nucleons and light charged 
particles Q36]], at variance with the partitions ob- 
served in the INDRA data which have a larger 
number of fragments. In this sense the "critical" 
values of [|5j] would rather be assimilated to "lim- 
iting" temperatures, or excitation energies, as de- 
scribed in [J37p. Above limiting values, a nucleus 
is unstable due to Coulomb forces, and can only 
de-excite through multifragmentation. 

What additional information is brought by the 
other signals quoted in table I? The universal 
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Table I: Summary of the phase transition signals described in the previous sections. Si„ c , eth, £rad, stand for the incident, 
thermal and expansion energies in AMeV (see text for precisions). SD means signal of spinodal decomposition. A 
indicates the exponent of the universal scaling, and the Bimod line gives the number of peaks in the probability distribution 
(© means two equal peaks). A "-" means that no measurement of the concerned observable was done. 



system 


Ni+Au 


Xe+Sn 


Xe+Sn 


Ni+Au 


Xe+Sn 


Xe+Sn 


Ni+Au 


Ni+Au 


£inc 


32. 


32. 


39. 


52. 


45. 


50. 


74. 


90. 




5.0 


5.0 


6.0 


6.5 


6.6 


7.0 


7.0 


7.5 


^■rad 


0. 


0.6 


1.2 


0. 


1.7 


2.0 


0. 


< 0.5 


SD 


no 


yes 


yes 


yes 


no 


no 






c<0 


yes 


yes 


yes 


yes 


no 


no 


yes 


yes 


A 




1/2 


1 




1 


1 






Bimod. 


© 


2 


© 


1 or 2 


2 


2 




1 



fluctuations are small for e t h < 5 AMeV, and 
large above (A = 1). This can be interpreted as 
the transition from the dominance of fusion evap- 
oration to that of fusion-multifragmentation pro- 
cesses around this energy; indeed at 32 AMeV 
the selected central collisions in Ni+Au present 
essentially the characteristics of an evaporation 
process [[24]]. In this sense, large fluctuations of 
Z max would sign the break-up of the system in- 
side the coexistence region. Bimodality of a pa- 
rameter related to the density of the system goes 
roughly in the same direction. Two peaks, indi- 
cating that the systems break in the coexistence 
region, are observed in all the cases except for 
Ni+Au at 90 AMeV This would indicate a pure 
gas phase, and contradict somewhat the measure- 
ment of a negative heat capacity. The equal abun- 
dance of the two phases occurs for e t h < 5-6 
AMeV, above that energy nuclear systems would 
mainly multifragment. 



Xe, Gd, Ta, U projectiles. 

Next, enriched knowledge of nuclear matter 
properties will come from taking into account 
its two component nature (neutrons and protons), 
by exploring collisions induced by exotic projec- 
tiles when they are available. Isospin fraction- 
ation should allow to distinguish the liquid and 
gas phases. 
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